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Abstract 

We consider the massless limit of Higgs gravity, where the graviton becomes massive when the 
scalar fields acquire expectation values. We determine the Vainshtein scale and prove that massive 
gravity smoothly goes to General Relativity below this scale. We find that the Vainshtein scale 
depends on the particular action of scalar fields used to give mass to the graviton. 
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I. INTRODUCTION 

In the recent paper [1] we (A.Ch.,V.M) have proposed a Higgs mechanism for gravity. 
In our model the graviton becomes massive as a result of spontaneous symmetry breaking, 
where four scalar fields acquire non- vanishing expectation values. As a result, three out of 
four degrees of freedom of scalar fields are absorbed producing a massive graviton with five 
degrees of freedom, while one degree of freedom remains strongly coupled. Our model is 
explicitly diffeomorphism invariant and, in distinction from bigravity theories, it is simply 
given by General Relativity supplemented with the action of four extra scalar fields. There- 
fore it is completely analogous to the standard Higgs mechanism used to give masses to the 
gauge fields, where masses are acquired as a result of the interaction with external classical 
scalar fields. For instance, in the standard electroweak theory one also uses four (real) scalar 
fields to give masses to three vector bosons, and one remaining degree of freedom becomes 
a Higgs boson. However, in distinction from electroweak theory, in our case the analogue of 
the Higgs boson remains strongly coupled and hence completely decouples from gravity and 
other matter. 

The theory with four scalar fields was exploited before by several authors (see 15Ml7| and 
references therein). In our case we have found the Lagrangian which resolved the problems 
that faced finding a consistent theory for massive gravitons. On one hand the model produces 
a graviton mass term with explicitly invariant form even for finite diffeomorphisms, and on 
the other hand, keeps the dangerous mode which could produce a ghost, in the strong 
coupling regime where it is completely harmless. In the linear order the mass term is of 
the Fierz-Pauli form which is uniquely fixed by the requirements of the absence of 
extra scalar degree of freedom. The analysis by Deser and Boulware js| however lead to 
the conclusion that in the massive theory the extra scalar degree of freedom reappears at 
nonlinear level and does not decouple, thus making massive gravity to be an ill-behaved 
theory. In distinction from [3], where diffeomorphism invariance is explicitly spoiled, our 
theory is diffeomorphism invariant and therefore the go a components of the metric remain 
always the Lagrange multipliers, while as we will show later, the scalar fields are always 
in the strong coupling regime above so called Vainshtein energy scale. This corresponds to 
extremely small energy and therefore the possible ghost is irrelevant. 

There were many interesting attempts to extend massive gravity beyond the linear ap- 
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proximation in a way where one can avoid the extra mode and ghost, also at the nonlinear 



level (see, for instance, \U 
recent interesting papers 
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24( and references there). In particular, in the 



24J | it was found an extension of the Fierz-Pauli action for 
which the ghosts are absent even at nonlinear level at the decoupling limit. 

The main purpose of this paper is to investigate the existence of a smooth limit of our 
model to Einstein gravity, when the mass of the graviton vanishes. It was noticed long ago 
by van Dam, Veltman and Zakharov j^, \^ that in linearized massive gravity the extra scalar 
mode of the graviton did not disappear and remained coupled to matter even in the limit of a 
vanishing graviton mass. In turn, this spoils predictions of General Relativity either for the 
perihelion precession or deflection of starlight. This effect is known as the van Dam-Veltman- 
Zakharov (vDVZ) discontinuity and was first thought to be a no-go theorem for massive 
theories of gravity . However, it was pointed out by Vainshtein that the discontinuity 
could be an artifact due to the breakdown of the perturbation theory of massive gravity in the 
massless limit jfj]. He has shown that in the case of gravitational field produced by a source 
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of mass Mo the nonlinear corrections become important at scales r < Ry = M m g 
(in Planck units) and conjectured that in the strong coupling regime General Relativity 
is restored. When the mass of the graviton m g vanishes the Vainshtein radius Ry grows 
and becomes infinite, thus providing a continuous limit to General Relativity in case the 
Vainshtein conjecture is correct. At distances r <C Ry, around a static spherically symmetric 
massive source of mass Mq the full non-linear strongly coupled massive gravity has to be 
considered in order to recover the Einstein theory, which makes the proof of the Vainshtein 
conjecture non trivial. The question of continuous matching of the solutions below and 
above the Vainshtein radius have been extensively addressed in recent literature. The first 
model where such a transition was demonstrated is Dvali-Gabadadze-Porrati (DGP) model 



which imitates many features of massive gravity 13|, |20j. There was a claim that in the 



bigravity version of massive graviton the corresponding solutions do not match but it 
was recently shown that this claim is not justified js-10|. 

In this paper we will find the Vainshtein scale and will prove Vainshtein conjecture in the 
Higgs model of massive gravity in the case when the gravitational field is produce by a source 
of mass M . Moreover, we will find how the concrete value of the Vainshtein scale depends on 
the nonlinear extension of the Pauli Fierz term, or in other words on the interactions of scalar 
fields used to produce massive gravity. As a result we will determine possible Vainshtein 
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scales for a wide class of Higgs gravity models. We will also derive in our model the leading 
corrections to the gravitational potential within Vainshtein scale, which are similar, but not 
identical to this type of correction obtained in the framework of the DGP model in 

Finally, we will discuss the implications of our results obtained in classical theory when 
extended to quantum theory. In particular we argue that in quantum theory there must be a 
cutoff scale at energies m g , above which the scalar fields enter strong coupling regime and 
completely decouple from gravity and other matter. Because this scale is extremely small 
for the realistic mass of the graviton it makes the problem of ghost which could appear 
only below this scale completely irrelevant. For the scalar and vector modes of the massive 
graviton the cutoff scale is an analog of the Planckian scale for the tensor graviton modes, 
which also become strongl 
agreement with results of 



y coup led above Planck scale. The obtained cutoff scale is in 

HQ- 
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II. HIGGS FOR GRAVITON: BASICS 



We employ four scalar fields (f> A , A = 0, 1, 2, 3 to play the role of Higgs fields. These will 
acquire a vacuum expectation value proportional to the space-time coordinates, thus giving 
mass to the graviton. Let us introduce the "composite metric" 



H 



AB 



(1) 



which is scalar with respect to diffeomorphism transformations. The field indices A,B, - ■ ■ , 
are raised and lowered with the Minkowski metric 7/ab- The diffeomorphism invariant action 
which will be used as our model, is given by 



S = -- d^x^gR^ — I d^xy/^g 




- v 2 H£H 



B 



where 



B 



(2) 



(3) 



is the traceless part of the "composite metric" and where we have set 87rG = 1. The param- 
eter v controls the symmetry breaking scale. As will be seen later, the induced mass of the 
graviton is equal to m g = Mv 2 and hence when v — > gravity becomes massless. It is clear 
that in this limit the only surviving term in action (T5]) is Einstein gravity and M 2 H A for 
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the four scalar fields, which are in the regime of strong coupling and do not possess linear 
propagators. In the phase with restored symmetry the total number of degrees of freedom 
is six: two of them describe massless graviton and four correspond to scalar fields which are 
decoupled from gravity at linear level. 

We show next that when the symmetry is broken, three out of four scalar fields are "eaten" 
and produce the massive graviton with five degrees of freedom, while the "surviving" degree 
of freedom will remain strongly coupled. In case when t)/0, the unique Minkowski vacuum 
solution of the equations of motion, g^ v = rj^, corresponds to the fields, which linearly grow 
with coordinates, that is, (fi A = y/vS^x 13 . Let us consider perturbations around Minkowski 
background, 

tf» = if + <p A = Vv~ (x A + X A ) (4) 

and define 

hi = \h$ -5 A = h A + d A X B + d BX A (5) 
+ d cX A d c XB + h A ,d c X B + h c B d cX A + h c D d D XB d c X A , 

where indices are moved with the Minkowski metric, in particular, xb = VbcX C an d h B = 
VBcfi A fiv W ■ We point out that we have included a factor y/v as coefficient of x A to obtain 
simpler expressions. In reality in all our results that will subsequently follow we have to 
make the replacement 

„ i /My , 



x A -> x A ^ =1 — 1 x' 

y/v \m g J 

This, however, will not effect most of our conclusions, and we will thus comment on it only 
when necessary. With the help of the expressions 



H = v (h + 4) , H a H B a = (h A B h B A - 



we can rewrite the action for the scalar fields in the following form 



MV 



Sa, = — - — / d xy/^g 



h 2 - h A h B A + ^¥ + ^h 4 



3 73 3 

— h -\ 



(6) 



We would like to stress that we did not use any approximations to derive (jHJ), and h B are 
diffeomorphism invariant combinations of the scalar fields and metric up to an arbitrary 
order. 
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III. PHYSICAL DEGREES OF FREEDOM OF THE MASSIVE GRAVITON 



We consider now small perturbations of the metric and scalar fields and neglect higher 
order terms. In this case 

h A B = h A B + d A XB + d BX A + C{h\ x 2 ), (7) 

and in the leading order, action (J6J) describes Fierz-Pauli massive gravity, where the mass 
of the graviton is equal to m g = Mv 2 . However, we have to stress that in distinction from 
the Fierz-Pauli theory our model does not break diffeomorphism invariance and coincides 
with this theory only in the unitary gauge where all x A = 0. In turn, imposing these gauge 
conditions completely fixes the coordinate system making the interpretation of the results 
rather obscure. If one would try to treat \ A as Stiickelberg "vector" field and consider 
the diffeomorphism transformations for the vectors rather than some obscure "fictitious" 
symmetries, then one unavoidably would conclude that the "vector components" must be 
treated as the perturbations of four scalar fields with nonzero background values, thus arriv- 
ing at our model. As we will see in the next section the difference between the noncovariant 
Fierz-Pauli approach and our model becomes even more dramatic at higher orders. However, 
we first study the linearized theory using Lorentz-violating approach to explicitly reveal the 
true physical degrees of freedom of the massive graviton. Namely, we use the method usually 
applied in cosmological perturbation theory and classify the metric perturbations according 



to the irreducible representations of the spatial rotation group [25|. The h 00 component of 
the metric behaves as a scalar under these rotations and hence 

h 00 = 20, (8) 

where is a 3-scalar. The space-time components h^ can be decomposed into a sum of the 
spatial gradient of some 3-scalar B and a vector Si with zero divergence: 

h 0i = B ti + S u (9) 

where = dBjdx 1 = diB and d z Si = 0. 
In a similar way hy can be written as 



hi, = 2ip5 i:j + 1E M + Fij + F jti + h i:j , (10) 
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where d l Fi = and d l hij — = h\. The irreducible tensor perturbations hij have two 
independent components and describe the graviton with two degrees of freedom in a dif- 
feomorphism invariant way. The scalar perturbations are characterized by the four scalar 
functions <p,ip,B, and E. In empty space they vanish and are induced entirely by matter, 
which in our case are the scalar fields. The vector perturbations of the metric Si and Fi are 
also due to the matter inhomogeneities The matter perturbations can also be decomposed 
into scalar and vector parts: 

X° = X°, X l = X l + ir fi (11) 
where diX 1 = 0. In the linear approximation, scalar, vector and tensor perturbations are 
decoupled and can be analyzed separately. 

Scalar perturbations. Up to first order in perturbations we have h a/3 = —Tf^rf^h^ 
and using the definition of hg in (JHJ) we find that in the leading order approximation 

( s )h° = -20 + 2x°, Wj,o = _ B . _ fr. + x o ? (5)^ = 2tfj6ik + 2E . k + 27r ^ (12) 

Substituting these expressions in (JSJ), keeping only second order terms, and expanding the 
Einstein action up to second order in metric perturbations we obtain the following action 
for the scalar perturbations: 

^S 2 S = J d 4 x |-3^ 2 + + [2A^ - m 2 (3^ + A(E + tt))] 

+ 2^A (b-e)+ m 2 [3^ (ip + X°) + (ty + X°) A(E + tt) 
+i(x -fi-7V) j ( X -fi- 7 r)J}, (13) 

where m 2 = M 2 v A and the dot denotes derivative with respect to time. We see that is a 
Lagrangian multiplier which implies the constraint 

2 

771 

= -y(3^ + A(£ + vr)). (14) 

Another constraint is obtained by variation with respect to B: 

iP = -^(x°-B-ir). (15) 

To simplify further the calculations we select the longitudinal gauge B = E = 0, which 
when used in conjunction with f fl4|) . simplifies the action (fTB"]) to 



-3t/> 2 + 



^S 2 S = J d 4 x 

+m 2 g (3V (V- + X°) + (2V> + *°) Atf + f ( X ° - tt) . ( X ° - tt) 
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Using constraints (JHJ) and (|T5|) with B = E = 0, imply 

rn^Avr = (2A - 3m 2 g ) if; (16) 

m*A X ° = - (2A + 3ml) 4 ( 17 ) 
which can be inverted to express n and x° m terms of if>: 

' - fc| - s) *• (18) 

Substituting these relations in the action above we obtain 



^S 2 S = / d 4 x 



6 • 4 

-3^ + + ( — if) 2 jif>,iif>,i - 3^ 2 

» m g rrig 



= -3 J d 4 x [if) (dt - A + m 2 g ) ip\ . (20) 

Note that the potential if) is gauge invariant with respect to infinitesimal diffeomorphism 
transformations: x a — > x a = x a + £ a . Therefore the derived result does not depend on the 
particular gauge we used to simplify the calculations of the action. First of all we see that 
the scalar mode which was non-propagating in the absence of the scalar fields has become 
dynamical. The variable u = \/6ip is the canonical quantization variable for the scalar degree 
of freedom of metric perturbations. It is entirely induced by perturbation of the scalar fields 
7r and x°- I n t ne linear approximation we have to be careful in taking the limit m g — > 
because of the inverse mass dependence in the relations ([IB]) and (fl9|) . In reality we have 
to consider instead equations f lTB]) and (fT7j) which implies that if) = as in the vacuum 
case. Thus the famous vDVZ discontinuity J, |5|] is not present. In addition, as mentioned 

before, when taking the limit m g — > we have to replace the fields n and x° with j 4 n 

1 

and (j^ J 4 x° but this leads to the same result that ■0 = 0. We note, however, that in the 
trig — > the Higgs action reduces to the M 2 H 4 term, and there are higher order non-linear 
contributions to if). In the next section we will show that above a certain energy scale the 
scalar mode ceases to propagate and becomes confined due to nonlinear corrections to the 
equations. As a result the vDVZ discontinuity is avoided completely and we obtain a smooth 
limit to General Relativity when symmetry is restored and the graviton becomes massless. 
Vector perturbations. For the vector perturbations 

(vy h o = _ s . _ ^ (v)£i =F(tk + Fjhi + jk + ^ ( 2i) 
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Up to second order in perturbations the action for the vector modes is 



iv) 6 2 S = j I d\v 



Fj, — Si) I Fi — Si 

k V / ,k 



+m 2 g ((^ + s^ (r + s t ) - + r), k + r), 



(22) 



Variation of this action with respect to Si gives the constraint equation 



A F - S, 



-m 



2 (r + s^ 



which allows us to express Si as 

Si 



A — Trig 



AFi + mli?) . 



(23) 



Substituting this expression into (|22|) we obtain 



m 2 „A 



[(Fi + r) (d 2 -A + m 2 g ) (F i+ f)]. 



(24) 



2 J 2 (A - raj) 

In the limit m 9 — )■ the action for the vector modes vanishes even after replacing x 1 — y 
— ) 4 x*. The canonical gauge invariant quantization variable in this case is the 3- vector 



V 1 



m 2 n A , .. 

{F t + r) 



(25) 



2(A-mJ) 

which describes two physical degrees of freedom as this vector satisfies an extra condition 
diV 1 = 0. 

Tensor perturbations. For the tensor perturbations the result is straightforward 



(T) 52 S = -t I d \ v 



kj {dt - A + m 2 g ) K 



(26) 



This action describes the pure gravitational degrees of freedom which have become massive. 
Because hij satisfies four extra conditions d l hij — = h\ the tensor perturbations have two 
physical degrees of freedom. 

Thus, we have decomposed the massive graviton with five degrees of freedom into physical 
gauge invariant components: a scalar part ip (with one degree of freedom), a vector part V % 
(2 degrees of freedom) and a tensor part hij (2 degrees of freedom). After quantization they 
acquire their independent gauge invariant propagators. 
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The metric components are the subject of minimal vacuum quantum fluctuations. In 
particular, the amplitude of the vacuum fluctuations of ip and hij at scales A <C l/n^g are 
about 

ip ~ hy ~ -, 

in Planck units. They become of the order of one at the Planck scale Zpi ~ 10 -33 cm where 
non perturbative quantum gravity becomes important. The amplitude of the vector vacuum 
metric fluctuations is much smaller. In fact, for A <C ^/ m g, their amplitude in the gauge 
Si = is scale independent and is equal to 

{v) h i:j ~ Fi tj ~ m. 

These results are valid only in linearized theory. While the result for the tensor fluctuations 
remains the same, we will show in what follows that the scalar and vector modes reach the 
strong coupling regime at the energy scale which is much below the Planck scale. 



IV. VAINSHTEIN SCALE AND CONTINUOUS LIMIT 



Let us first consider how the static interaction between two massive bodies is modified 
in the Higgs model with massive graviton. In quantum field theory this interaction is 
interpreted as due to the exchange by gravitons with corresponding quantum propagators. 
This interpretation is very obscure from the physical point of view because the Newtonian 
force is not directly related to the propagation of gravitons. It is, however, the price to be 
paid in order to preserve explicit Lorentz invariance of the theory. In our approach one does 
not need to go to quantum theory to answer this question. The interaction is entirely due 
to the static potentials and ip which are present due to the massive body. Let us take the 



Newtonian gauge 25], where B = E = so that the metric takes the form 

ds 2 = (1 + 2(f)) dt 2 -(l- 2^0 5 ik dx i dx k (27) 

First we have to derive the equations that this metric should satisfy in massive gravity. We 
consider only static solutions so all time derivatives vanish and action (fl3l) simplifies to 



^5 2 S 



J d 4 x {ifijifjj + [2A^j - m? g {3i/; + An) - T 00 ] 



+m l[3^ + 2ijAn + \x l X°i]}- (28) 
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We have added a term which describes the interaction with an external source of matter for 
which only the T 00 component of the energy momentum tensor does not vanish. Varying 
this action with respect to 0, -0, X° an d tt we arrive to the following equations: 

A^ = ^(3^ + Att) + — , Aty-<i>-m 2 g v) = 0, (29) 

A X ° = 0, A(2^-0) = O. (30) 
It immediately follows from fl30|) that x° — and ip = 0/2, while equations (129]) simplify to 

A(0 + ^) = 3m^ + T 00 , (31) 

or taking into account that ip = 0/2 we obtain 

4 /T 00 \ 

(A-^)0 = ^(V)- (32) 
For the central source of mass Mo the solution of this equation is 

4 M n 4 

<P = ~\— e- m ° r = -<P N e- m ° r , (33) 
3 r 3 

where <Pn = — Mo/r is the Newtonian gravitational potential. At scales r <C 1/^g the metric 
takes the form 

ds 2 = (l + | (20^ rft 2 - U - ^>jv) 5 ik dx i dx k . (34) 

The bending of light is determined by the + -0 combination of the metric components. In 
General Relativity, where -0 = 0^, this combination is equal to 20^. In the case of massive 
gravity 

4 2 

4> + i> = ^<PN + g0iV = 20AT, (35) 

i.e. we obtain the same prediction for the bending of light. However, the gravitational 
potential which, for instance, determines the motion of planets has increased by factor 4/3 
compared to the Newtonian potential, independently of the mass of the graviton. This extra 
contribution survives even in the limit of zero mass. If one would redefine the gravitational 
constant to get the correct Newtonian potential then obviously the bending of light would 
be wrong. This is a manifestation of vDVZ discontinuity, which in quantum field theory is 
interpreted as due to the propagation of the extra scalar mode in addition to the two tensor 
degrees of freedom. Because this scalar mode is coupled to the trace of the matter the result 
remains unchanged for photons, but changes by the corresponding factor for non-relativistic 
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matter. Note that we have re-derived this result in a purely classical theory without any 
reference to the tensor degrees of freedom or the "true" graviton. 

The paradox with vDVZ discontinuity, which implies that the graviton must be strictly 
massless was resolved when Vainshtein found a new scale Ry in massive gravity and sug- 
gested that for r < Ry the scalar mode decouples and General Relativity is restored. 

We will now show how this happens in our theory, and prove that General Relativity is 
smoothly restored below the Vainshtein scale. For that we will need to consider the higher 
order corrections to the action (1281) . First of all we notice that because the gravitational 
potentials with which we are dealing are always much smaller than unity, we can safely 
ignore the terms of order 3 ,0-0 2 etc. compared to <f) 2 , ... because they cannot change the 
solutions of the equations drastically. We will also ignore the terms <fi 2 (An) compared to 
<p (A7r) etc. because they are subdominant. Therefore, the only contribution to the higher 
order corrections which we will take into account will come purely from the matter scalar 
fields. In addition we will skip all terms with x° since they vanish in the leading order. 
Hence, the only relevant terms of the third order, which should be added to the action (1281) 
are: 



d x 



- (Ann, ik n^ ik - ■K ;ki 7r ii j7r j j k ) + — (An) 3 - - 
9 



+ 2ip) ir ik ir ik 



+2^ (Avr) 2 + — (3tp - 0) (Avr) 2 + O (tfj 3 , ifj 2 ^ V 2 Avr, (ptpAn...) 
16 



(36) 



These third order corrections modify the equations obtained by variation with respect to ip 
and n in the following way: 



A (ip - <p - m 2 g n) + 



m. 



3 , « 1 59 , . , 2 

2 VP ~ 2 V0 + 2 7r .** 7r .<* ~ 54 ( A?r ) 



0. 



and 



(37) 



A (2ip -<p) + (Atttt 
9 



ikj 



,ik + 2^ ( 7T ' ifc7r > ifc ^ 



ik 



+ —A (Any + O (<j), lk 7r, ikj A^Avr, ...) 
16 



0. 



(38) 



Equation (|38p is the main equation where non-linearities begin to play an important role 
allowing us to avoid the condition A (2ip — 0) =0, and thus resolve the problem of vDVZ 
discontinuity. In fact, this condition means that the scalar perturbations of the curvature 
must vanish, 5R = 0, and this was the main obstacle leading to the troubles with restoring 
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General Relativity in the limit of vanishing graviton mass in the paper j3j. Assuming that 
7Tjfc, An <C 1 (this assumption will be checked a posteriori), and keeping only the leading 
terms in equations ( 157)) and ( 155)) we obtain 

A (ij> - <f) - m 2 g ix) = 0, A (2if> - 0) + <9 6 vr 2 = 0, (39) 

where by d 6 ir 2 we denoted all quadratic 7r terms in ( 138)) . Using the first equation in (139)) to 
solve for A<p, the second one simplifies to 

A (ip + m 2 vr) + d 6 7i 2 = 0. (40) 

Taking into account that A ~ d 2 and estimating <9 6 7r 2 in spherically symmetric field as 
O (1) 7r 2 /r 6 , this equation becomes 

ip + m 2 g n + O (I) r-\ 2 ~ 0, (41) 

The behavior of 7r as a function of r crucially depends on whether the second or third term in 
this equation is dominating. To estimate the scale when both terms are comparable, which 
is called the Vainshtein scale Rv, we set 

m 2 7i ~ O (1) r _4 7r 2 ~ ip, 

and from here find that 



- Mr=R v = m t R V- 



(42) 



In the case of a gravitational field produced by the mass Mo in the vacuum ij) ~ — Mo/r, 
the Vainshtein scale is equal to 



R 



v 



Mn^ 1/5 



rrvt 



For r ^> Ry the last term in (1411) is small compared to the second one and we obtain 

5 N 



7T 



•1 + 



r 



(43) 



(44) 



In this limit the quadratic terms in the second equation in (139)) are negligible and from the 
first equation in (139)) we find that 

5 N 



ifj — (j) = —Ip 



1 - o 



Rv 

r 



(45) 
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This implies that in the leading order ip = 0/2 in complete agreement with the result which 
we have obtained above in linearized massive gravity. It is easy to check that the condition 
d 2 n <C 1 which we have used to simplify equations (157)) and (1551) is also satisfied. In fact, 

d\ ~ -4- ~ iv ~ (^Y7 -f-") 2 , (46) 



9- 



and hence <9 2 7r <C 1 for all r > Ry if i?y <l/m 

At scales smaller than Vainshtein radius, that is for r <C Rv the third term in (1411) is 
larger than the second one and hence 

2 2 

l + O 1 ^= + ... 



7r ~ (1) r v— if) 
= 0(1)-* ff 



- r / r \ 5 / 2 



(47) 



^-0 = 0(1)^(^1 +... (48) 



Using this expression in the first equation of (159"]) we then find that in the leading order 

5/2 

Ry, 

For r <C we find that if) = up to corrections of order if) (r/i?y) 5 ^ 2 . Because of d 2 n ~ 
\J—i\) the condition <9 2 7r <C 1 is always satisfied. The dominating quadratic corrections to 
equation (15T1) is of order m 2 ip ~ m 2 (<9 2 7r) 2 so they change only the mass term which is 
irrelevant within the Vainshtein scale. Taking into account that if> = for r <C i?y and 
neglecting the mass term, equation (15T1) in the leading order is reduced to 

A0 = — , (49) 
and thus General Relativity is restored within Vainshtein scales up to the corrections 

\ 5/2 

*;) • < 50 > 

[ — I 

which are much smaller than the corresponding corrections in DGP model [22]. One could 
ask whether any higher order corrections would be able to spoil the obtained results? The 
most dangerous of these corrections in every next order will come as the corrections to the 
previous order multiplied by d 2 7r <C 1. Therefore they are completely negligible. 

We now consider the implementation of our results derived classically in quantum field 
theory. In the explicitly Lorentz invariant approach the change of the interaction strength 
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at scales exceeding the Vainshtein scale is interpreted as due to exchange by the scalar mode 
if) of the massive graviton in addition to the two tensor modes of the massless graviton. As 
we will argue, this scalar mode becomes strongly coupled below Vainshtein scale and as a 
result completely decouples from the gravity and matter entering the confinement regime. 
This is similar to QCD, where the "soft" modes do not participate in the interactions of 
highly energetic quarks below the confinement scale. Although for quantum fluctuations 
one cannot neglect the time derivatives as in the static case, we can, however, estimate 
the time derivatives to be of the same order of magnitude as spatial derivatives and use 
the formulae derived for the static case. Keeping in mind that the amplitude of the scalar 
quantum fluctuations at the length scale A is about if) ~ 1/A from (112]) we obtain that at 
scales smaller than 

A s ~ m g " 4/5 , 

these scalar modes should be in the strong coupling regime, where nonlinear corrections 
cannot be neglected. Note that the metric fluctuations which are of order ip ~ m g still 
remain small at this scale. In distinction from the case when gravitational field is produced 
by an external source the estimate d 2 n ~ \J — ip is not justified for quantum fluctuations 
for A < A s . However, assuming that at the scales which are just a bit smaller than A s 
one can still use this estimate to find that the last term in action (1131) . which is of order 
<9 4 7r 2 ~ if), becomes dominant compared to the terms of order tp 3 ^ 2 and ip 2 . As a result the 
scalar mode if) loses its linear propagator and decouples, entering the strong coupling regime 
where nonlinear corrections will prevent its unbounded growth for every A < A s as m g — > 0. 
As a result the terms proportional to m 2 in the action ( Tl3|) will vanish and General Relativity 
is smoothly restored in this limit. A similar thing happens with the vector modes. Therefore 
in the limit m g — > only the tensor modes hit with two degrees of freedom survive. They 
enter the strong coupling regime at the Planckian scale. The energy scale A" 1 should be 
taken as a cutoff scale for the scalar mode if) of graviton in all diagrams where this scalar 
mode participates. Above this scale our scalar fields tt and x° which were producing the 
extra degrees of freedom for the massive graviton are also in the confined regime and the 
symmetry is restored. These strongly coupled fields are completely decoupled from gravity 
and the rest of the matter. In the case when the mass of the graviton is of the order of 
present Hubble scale the cutoff scale is extremely small of order 10 -18 eV. At higher energies 
the ghost, even if it would exist, completely decouples. Therefore the question about ghosts 
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at the nonlinear level becomes irrelevant. 



V. HOW UNIVERSAL IS THE VAINSHTEIN SCALE? 

The expression (H3j) for the Vainshtein scale was derived first in the case of Fierz-Pauli 
mass term which is unique in four dimensions, because only in this case there are no ghosts 
propagating at the linear level. We have obtained the same result in our Higgs model with 
the action fl2]). It is natural to ask whether it is the unique universal scale for all models with 
Fierz-Pauli mass term or it depends on a particular nonlinear extension of this term. Let us 
show that in our theory the Vainshtein scale, in fact, depends on the nonlinear completion 
of the theory and determine all possible extensions of the model which lead to different 
Vainshtein scales. With this purpose we first consider instead of ([2]) the following action for 
the scalar fields 

2 / TT \ 3 



S, = ^jd^^-9 "(f-l) + 43 / ? (f- 1 ) -H A B H B A +0{{H-Af) 



(51) 



where without loss of generality we have set the parameter of the symmetry breaking to 
unity. The terms O ((if — 4) 4 ) must be taken in such a way as to avoid the appearance of 
other vacua, besides if = 4. One can easily verify that there are infinitely many extensions 
of the required type. This action, when rewritten in terms of variables defined in (JSJ), 
with v = 1, takes the form 

S * = ~8 / ^ X ^ ^ ~ h ^ + ffi + (h 4 )} , (52) 

where m 2 = M 2 . For ^> 1, the main contribution to the cubic action fl36|) is of order 
j3 (A7r) 3 and the second equation in (|39|) is modified to 

A (2^> - 0) + 3/3A (Att) 2 = 0. (53) 
Then using the first equation in (|39l) and considering the spherically symmetric case we find 

ip + m 2 vr + O (1) I3r- 4 n 2 ~ 0, (54) 
and correspondingly the Vainshtein scale in this case is 

Rv ^{^4) ■ (55) 
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Thus, we see that taking large enough (3 in action (I5ip we can obtain an arbitrarily large 
Vainshtein scale for given masses of the source Mq and the graviton m g . 

Next we would like to address the question whether one can obtain a smaller Vainshtein 
scale compared to fl4"3"j) . For that let us first consider the action 



Sd 



M 2 



d x 



-g 



-6 



H 

T 



1 



H 

T 



\h a b h b a + 



+\h^hEh c a - \hh a h b a +0((H- 4) 4 ) 



(56) 



where the terms O ((H — 4) 4 ) are taken in such a way as to avoid the vacuum at H = 12. 
Rewritten in terms of h B , action fl56l) becomes 

,2 



m g 
Sa> = — — 



d x 



-fj 



¥ - h A B h B A + ~ (h A B h B h c A - h A B h B A h) + O (h*) 



(57) 



where m 2 g = M 2 . It is clear that in the lowest order it reproduces the Fierz-Pauli term, but 
in higher orders it is quite different from ((6]). The action (1571) concides with the action first 



derived in 
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24j from the requirement of the absense of ghost in decoupling regime up to 
the third order. If we consider the case of the static gravitational field we find that in the 
third order the action does not contain terms of the form <9 6 7r 3 . Hence, by keeping only the 
leading terms we find that the second equation in fl39l) will be modified to 

A (2ip -<f)) + dV = 0. (58) 

Considering the spherically symmetric case and using the first equation in (|39|) . which is 
still valid up to the leading order, we find that equation (1411) has to be replaced by 



V> + m> + 0(l)r- 6 7T 3 



0. 



(59) 



The Vainshtein scale will be determined by the condition that all three terms in this equation 
become comparable, that is, 

i> ~ m 2 g 7i ~ r" 6 7r 3 (60) 
and hence the expression determining this scale is 



m 3 g R 3 v . 



In particular, in the case of static field produced by mass Mq, we have 



R 



Mo 
m 3 



1/4 



(61) 



(62) 
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To obtain the correction to the Newtonian potential at r Ry we note that at these scales 

2, 



rr .-n. r 2 ^ 1 / 3 a nd use of the first equation in fl39|) leads to 

\ 8 / 3 



p , ^ (63) 

If we set the mass of the source in (1621) to be equal to the Planck mass, the corresponding 
cutoff scale in quantum theory for the decoupling of the scalar mode is obtained : A s = 

-3/4 

m g . 

In principle, there are enough different combinations of which can be added to the 
action (l5"Tj) to remove all the terms of the form (d 2 n) k for all k < n, so that the first survived 
terms of this structure are (<9 2 7r)™' . Notice that such action is unique up to the order h n . In 
this case, the Vainshtein scale is determined by the condition 

2Q-1) 

- *P\ r=Rv = (m g R v ) n-2 . (64) 
In the case of static gravitational field due to a massive source M this yields 

R v = (M^- 2 mf l - n) )^ (65) 
and the correction to the gravitational potential for r <C Rv is of order 

(66) 



r 

Rv 



in agreement with 



26| . In the limit when n — > oo the Vainshtein scale is Rv = M^ 3 m^ 2 ^ 3 . 
It coincides with the corresponding scale in the DGP model. However, the corrections to 
the gravitational potential which decay as (r/Rv) 3 seem different. In this limit the theory 
is unambiguous, but one could write it only as an infinite series. In turn this indicates that 
such theory is most probably nonlocal. Moreover, because d 2 ir 4 1 we completely lose 
control of higher order corrections and hence the results become completely unreliable. 



VI. CONCLUSIONS 



We have addressed the most fundamental question of all theories of massive gravity - 
can massive gravity be a consistent theory not contradicting to current experimental and 
theoretical knowledge? In this paper we have treated gravity mostly as a classical field 
theory and have explicitly investigated the issue of a smooth limit of massive gravity to 
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General Relativity. With this purpose we first determined the physical degrees of freedom 
of the massive graviton generated via Higgs mechanism. This was done in the framework of 
irreducible representations of the three dimensional rotation group, where the five degrees 
of freedom of the graviton are described in terms of a tensor mode with two degrees of 
freedom and vector and scalar perturbations due to the scalar fields. The propagator for 
each of these five constituents of massive gravity was derived separately. In the linear 
approximation the origin of the well-known vDVZ discontinuity at the zero mass limit was 
traced to the constraint equations and it was shown how the scalar and vector modes of 
metric perturbations become non-dynamical in this limit. 

It has been suggested long ago that the linear perturbation theory of massive gravity fails 
at length scales below the Vainshtein scale and one has to consider the full nonlinear theory 
to recover General Relativity below this scale. We have determined the Vainshtein scale in 
Higgs gravity, with Fierz-Pauli mass term, and found the explicit solution for the spherically 
symmetric gravitational field. We have shown that the massive gravity solution outside 
the Vainshtein scale smoothly goes to the General Relativity solution in the region deep 
inside the Vainshtein scale. Thus the classical results and predictions of General Relativity 
are recovered inside the Vainshtein scale and at distances exceeding the Vainshtein radius, 
massive gravity strongly differs from Einstein theory. This means that the scalar mode of 
massive graviton decouples at Vainshtein scale and enters the strong coupling regime. In the 
limit of vanishing mass, when Vainshtein radius becomes infinite, the symmetry is restored 
and our theory is reduced to General Relativity with four scalar fields which are confined and 
thus decoupled from gravity and other matter. Based on these results we have argued that 
in quantum theory there is a cutoff energy scale above which the scalar fields responsible for 
the scalar and vector modes of the massive graviton are strongly coupled and confined and 
hence harmless. For the realistic graviton mass this scale is extremely low. Therefore, the 
question about extra scalar mode and ghost instability seems to be irrelevant in our model. 

We have found how the Vainshtein scale depends on the particular Higgs model or, in 
other words, on the nonlinear extension of the Fierz-Pauli mass term. In particular, we 
have shown that for given masses of the graviton and source, the Vainshtein length scale 
depends on the Lagrangian of the scalar fields and can be made arbitrary large. On the other 
hand, we have also constructed Lagrangians, which produce smaller scales compared to the 
standard one. However, the smallest possible scale seems to be larger than M^rrig 2 ^ 3 . 
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Finally, we have calculated the corrections to General Relativity within the Vainshtein 
scale which could, in principle, be interesting from experimental point of view. 
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